Understanding the impact of ammonium ion substitutions on heterogeneous ice nucleation.

Understanding the mechanisms underpinning heterogeneous ice nucleation in the presence of ionic inclusions is important for fields such as cryopreservation and for improved models of climate and weather prediction. Feldspar and ammonium are both present in significant quantities in the atmosphere, and experimental evidence has shown that feldspar can nucleate ice from ammonium-containing solutions at temperatures warmer than water alone. In recent work, Whale hypothesised that this increase in nucleation temperature is due to an increase in configurational entropy when an ammonium ion is included in the ice hydrogen bond network (T. F. Whale, J. Chem. Phys., 2022, 156, 144503). In this work, we investigate the impact of the inclusion of an ammonium ion on the hydrogen bond network by direct enumeration of the number of structures found using Rick's algorithm. We also determine the energy of these systems and thus compare the effects of enthalpy and entropy to test Whale's hypothesis. We find that the inclusion of an ammonium ion increases the total number of configurations under conditions consistent with a realistic surface charge. We also find that the enthalpic contribution is dominant in determining the location of the ammonium ion within the structure, although we note that this neglects other practicalities of ice nucleation.

Interplay of multiple clusters and initial interface positioning for forward flux sampling simulations of crystal nucleation.

Forward flux sampling (FFS) is a path sampling technique widely used in computer simulations of crystal nucleation from the melt. In such studies, the order parameter underpinning the progress of the FFS algorithm is often the size of the largest crystalline nucleus. In this work, we investigate the effects of two computational aspects of FFS simulations, using the prototypical Lennard-Jones liquid as our computational test bed. First, we quantify the impact of the positioning of the liquid basin and first interface in the space of the order parameter. In particular, we demonstrate that these choices are key to ensuring the consistency of the FFS results. Second, we focus on the frequently encountered scenario where the population of crystalline nuclei is such that there are multiple clusters of size comparable to the largest one. We demonstrate the contribution of clusters other than the largest cluster to the initial flux; however, we also show that they can be safely ignored for the purposes of converging a full FFS calculation. We also investigate the impact of different clusters merging, a process that appears to be facilitated by substantial spatial correlations-at least at the supercooling considered here. Importantly, all of our results have been obtained as a function of system size, thus contributing to the ongoing discussion on the impact of finite size effects on simulations of crystal nucleation. Overall, this work either provides or justifies several practical guidelines for performing FFS simulations that can also be applied to more complex and/or computationally expensive models.

The seven deadly sins: When computing crystal nucleation rates, the devil is in the details.

The formation of crystals has proven to be one of the most challenging phase transformations to quantitatively model-let alone to actually understand-be it by means of the latest experimental technique or the full arsenal of enhanced sampling approaches at our disposal. One of the most crucial quantities involved with the crystallization process is the nucleation rate, a single elusive number that is supposed to quantify the average probability for a nucleus of critical size to occur within a certain volume and time span. A substantial amount of effort has been devoted to attempt a connection between the crystal nucleation rates computed by means of atomistic simulations and their experimentally measured counterparts. Sadly, this endeavor almost invariably fails to some extent, with the venerable classical nucleation theory typically blamed as the main culprit. Here, we review some of the recent advances in the field, focusing on a number of perhaps more subtle details that are sometimes overlooked when computing nucleation rates. We believe it is important for the community to be aware of the full impact of aspects, such as finite size effects and slow dynamics, that often introduce inconspicuous and yet non-negligible sources of uncertainty into our simulations. In fact, it is key to obtain robust and reproducible trends to be leveraged so as to shed new light on the kinetics of a process, that of crystal nucleation, which is involved into countless practical applications, from the formulation of pharmaceutical drugs to the manufacturing of nano-electronic devices.